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Aims. We searched for very high energy (VHE) y-ray emission from the supernova remnant Cassiopeia A 

Methods. The shell-type supernova remnant Cassiopeia A was observed with the 17 meter MAGIC telescope between July 2006 and January 2007 
for a total time of 47 hours. 

C Results. The source was detected above an energy of 250 GeV with a significance of 5.2 cr and a photon flux above 1 TeV of (7.3+0.7 j,a,±2.2„,.i) x 
lO^'^cm^-i-"' . The photon spectrum is compatible with a power law dN/dE ocE"'^ with a photon index r=2.3+0.2j,a,±0.2svj. The source is point-like 
ILJ within the angular resolution of the telescope. 
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ABSTRACT 
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1. Introduction (Il997l) . The progenitor's init ial mass was large, estimated to be 

between 15 and 25 Mq, see ' Young et al.l ([2006). The morphol- 

Cassiopeia A (Cas A), with right ascension (RA) and declination ^^e remnant as seen in optical. X-ray and IR wavelength 

(DEC) (23.385 ',58.800°), is a prominent shell type supernova consists on a patchy and irregular shell with a diameter of 4' 

remnant and a bright source of synchrotron radiatiori observed at (4 at 3.4 kpc). The supernova blast wave is expanding into a 

radio frequencies, see Bell et al. ('1975'); Tuffs et al. ('1986'), and ^^^^ bubble formed from the previous wind phases of the pro- 

in the X-ray band, see Allen et al. ( 1997), Favata et al. ( 1997). ^^^^^^j. star; this plays an impo rtant role in shock acceleration of 

The remnant results from the youngest known Galactic super- see lBerezhko et all (l20()3b 
nova, whose explosion t ook place around 1 680. Its distance was 

estimated at 3.4 kpc bv lReed et al.l (Il995l). High resolution X - At TeV energies, Cas A was detected by the HEGRA 

rav images from the Chandra satellite, see 'Hu ghes et al.l(l2000b . Stereoscopic Cherenkov Telescope System, which accumu- 

reveal a shell-type nature of the remnant and the existence of a lated 232 hours of data from 1997 to 1999. TeV y- 

central object. The progeni tor of Cas A. was prob ably a Wolf- ray emission was detected at 5cr level and a flux of 

Rayet star, as discussed in iFesen et al.l (Il99ll) and llvudin et al.l (5.8 + l.2stat± \.2svs) 10^'"^ phcm~^s~' above 1 TeV was derived, 

as discussed in lAharonian et al.l (12001 1) . The spectral distribution 

Send offprint requests to: between 1 and 10 TeV was found to be consistent with a power 

E. Ofia-Wilhelmi, e-mail: emmaOapc. univ-paris7.fr law with a differential spectral index of -2.5+0.4.„fl,+0.1.5 
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per limits were consistent with the HEGRA detected flux level. 
At lower energy, EGRET set an upper limit for a flux below 
12.4 X 10"** cm'^s'\ see Esposito et al. ( 1996). 

The HEGRA detection makes Cas A a good scenario to test 
the supernova remnant emission at lower energies, in particular 
for trying to further distinguish between leptonic and hadronic 
models for the origin of the y-ray emission. A summary of the 
observations and analysis results is given in Section|2] the results 
are reported in Section [3] and finally a comparison of MAGIC 
detection with the existing model predictions for the TeV y-ray 
emission on Cas A is discussed in Section]?] 

2. Observations and Data Analysis 

The MAGIC (Major Atmospheric Gamma Imaging Cherenkov) 
Telescope is located on the Canary Island La Palma (2200 m 
asl, 28°45W, 1 T54'W) and has a 17 m-diameter tessellated re- 
flector dish, see Cortina et al. ] (120051) . The total field of view is 
3.5°. The accessible energy range spans from 50-60 GeV (trigger 
threshold at low zenith angle) up to tens of TeV. The telescope 
angular resolution (sigma of the Gaussian fit to the point spread 
function, PSF, o-psf) is about 0.09°. 

Cas A observations were performed between June 2006 and 
January 2007 for a total observation time of 47 hours after qual- 
ity cuts, namely, after rejecting runs with detector problems or 
adverse atmospheric conditions. The zenith angle ranged from 
29° to 45° and averaged 35°. The observ ation technique ap- 
plied was the so-called wobble mode, see iDaum et all (Il997h 
in which the telescope pointed alternatively for 20 minutes to 
two opposite sky positions at 0.4° off the source. Most of the 
data were taken under moderate moonlight illumination (86% 
of the scheduled observation time). Depending on the different 
moonlight levels, the resulting PMT anode currents ranged be- 
tween 1 jjA and 6 /uA, as compared to a typical anode current 
of 1 juA for dark night observations. Correspondingly, the trig- 
ger discriminator threshold (DT) was varied between 15 and 30 
arbitrary units (a.u.) to keep a low rate of accidental events. 
The mean trigger discriminator threshold during the observa- 
tions was 19 a.u, which corresponds to 13.3 photoelectrons (PE). 
Briefly, the impact of the rise of DT can be summarized by 
a decrease on the relative y-ray efficiency from 1 (dark obser- 
vations) to 0.84 while the relative sensitivit}[[] worsens from 
2.5 to 2.7% with respect to the Crab flux. Although this ef- 
fect is important for images containing a low number of PE's 
(low size), the energy threshold rise ( ~ 5 GeV) is negligible 
compared to the rise due to the medium to high zenith angle. 
Hence, the moderate moonlight illumination did not substan- 
tially reduce the telescope perform ance. More details on the 
Moon data analysis are discussed bv Albert et alj (l2007al) . Dark 
and Moon data were analyzed together using the standard anal- 
ysis and calibration programs for the MAGIC collaboration, e.g. 
iGaug et al. (2005). The images were cleaned using absolute tail 
and boundary cuts of 10 and 5 PE, respectively. For the y/hadron 
shower separation, the shower images w ere pa rameterized using 
the Hillas parameters, see Hillas et aTj (Il985h . These variables 
were combined for y/hadron separation b y mea n s of a Random 
Forest classification algorithm, see .Bock et al.l (l2005l) . trained 
with MC simulated y-ray events and data from galactic areas 
near the source under study but containing no y-ray sources. The 
Random Forest method calculates for every event a parameter 
dubbed HADRONNESS (H), which parameterizes the purity of 

' Minimal flux detectable with 5cr significance in 50 hours of obser- 
vations. 
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Fig. 1: The upper panel shows the distributions of 6- (measured in 
degree^) in the direction of the source (black dots) and anti-sources 
(blue shaded histogram). A lower SIZE cut of 400 PE was applied. The 
lower panel shows the 0^ distribution after background subtraction. The 
vertical line shows the optimum angular cut. The red function corre- 
sponds to the telescope PSF and is derived from Crab nebula data taken 
in the same observational conditions as those for Cas A. The black dis- 
tribution is the result of a Gaussian fit to the excess distribution. 



hadron-initiated images in the multi-dimensional space defined 
by the Hillas variables. 

The 9^ distribution is computed for the source position, 
where 9 is the angular distance between the source position in 
the sky and the reconstructed origin position of the shower. The 
reconstruction of individual y-ray arrival directions makes use of 
the so-called DISP method ( Domingo-Santamaria et al.1 (120051) V 
The expected number of background events are calculated using 
five regions symmetrically distributed for each wobble position 
with respect to the center of the camera and refered to as anti- 
sources. 

The optimum H and the angular cuts were derived using dark 
night Crab data of the same epoch and in the same observation 
conditions (zenith range, astronomical nights). The use of a dark 
night data sample in optimizing the telesc ope sensitivity is justi- 
fied by the results in lAlbert et al.l (l2007al) . For the spectral anal- 
ysis, the energy of each individual y-ray candidate was also es- 
timated using the Random Forest technique. The average energy 
resolution for the analyzed energy range was 20%. 

3. Source detection, extension and energy 
spectrum 

The so-called 9^ distributions for the source and anti-source po- 
sitions are shown in Figure JT] for a lower SIZE cut of 400 PE, 
which optimizes the MAGIC signal to noise ratio. The black 
points correspond to the source position whereas the blue-shaded 
histogram coiTesponds to the anti-sources. The subtraction of the 
two histograms shows the excess in the direction of Cas A. An 
excess of Ne.uess-i5 7 with a signifi cance of 5.2 cr (using the 
likelihood method ofl Li& Mai (119831) ) is detected within the re- 
gion 0. 13° centered at the HEGRA position. 

Figure ]2] shows the excess map of y-ray candidates with 
images larger than 400 PE. The map has been smeared with a 
Gaussian of cr=0.07''. The source position has been determined 
by ways of a fit of the non-smeared sky map to a bidimensional 
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Fig. 2: Sky map around the position of Cas A. A lower cut in SIZE of 
400 PE was applied. The green crosses mark the 2 wobble positions. 
The red cross indicates the MAGIC best fit position. The black cross 
marks the HEGRA source position, which is within 1 standard deviation 
from the MAGIC one. The bars of the crosses for both the MAGIC and 
HEGRA marks correspond to 1 sigma statistical errors. 



Gaussian function. The best fit position coordinates are RA - 
23.386+0.003 j,„,+0.001.„,, h and DEC = 58.81+0.03,,„,+0.02°,, 
(for more details on the s ystematic uncertain ties in the source 
position determination, see lBretz et a n (l2005h V 

In X-rays and radio-frequencies Cas A has an angular diame- 
ter of 0.08°, which is just on the limit of the MAGIC angular res- 
olution. The MAGIC system PSF is derived from MC simulation 
for a point source, and is found to be cr^„y =0.090±0.002° (shown 
in Figure |2|. Thi s value was validate d with Mkn 421 and Crab 
Nebula data (see lAlbert et alJ (l2007bh ). To further constrain the 
extension of the source we fit the excess with a Gaussian func- 
tion convolved with the PSF (F-A ■ exp{-0.5 9^l(_cr]^c + "^^i/)))- 
We obtain a value of (Jsrc which is compatible with zero within 
the fit error Figure [T] shows the telescope PSF and the result of 
the Gaussian fit (dotted blue curve). 

Figure [3] shows the reconstructed spectrum above 250 GeV. 
The spectrum is consistent with a power law (dN/dE oc E""^). 
The diff^erential flux at 1 TeV is (1.0+0.1j,n,+0.3„,s) x 10"'^ 



TeV" cm" s with a photon index of r=2.4±0.2j,a,+0.25,.5.The 
systematic eiTor is estimated to be 35% i n the flux level dete r- 
mination and 0.2 in the spectral index (see lAlbert et al.l (l2006l) ). 
The measur ed spectrum was u nfolded using the Gauss-Newton 
method, see ISchmelliiigl ( 11994 . The ;t-^/d.o.f of the fit is 2.83/3. 
The 1 cr error limit on the flux fitted is also added as a grey band. 
The Cas A flux corresponds to an integral flux above 1 TeV of 
3% of the Crab nebula flux above the s a me thr eshold (in red 
dashed line in figure [3] see [Albert et all (l2007bl) ). The Cas A 
spectrum measured by HEGRA is also shown as a blue solid 
line. The spectrum measured about 8 years later by MAGIC is 
consistent with that measured by HEGRA for the energies above 
1 TeV, i.e, where they overlap. 



4. Discussion 

The VHE MAGIC 47-hour observation of Cas A confirms the 
source detection by HEGRA after a multi-year integration of 
232 hours and at the same time significantly extends the en- 
ergy spectrum down to about 250 GeV. Cas A is detected with 
more than 5 cr at a flux level compatible with the HEGRA mea- 
surement for those energies explored in common. The difer- 
ential flux at 1 TeV measured by MAGIC is 1.0+0. ls,„, x 
2Q-i2y-gy-i^^-i^-i compared with the one measured by 
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Fig. 3: Cas A spectrum above 250 GeV. The blue line represents the ear- 
lier measurement by HEGRA. The red line represents the Crab nebula 
spectrum. The shaded area is the Icr statistical error of the fit. 



tween the two measures is excellent not only in the determina- 
tion of the flux level but also in the spectral index measured. 
Although the eiTors in the spectral index are large, there is no ev- 
idence for a high energy cutoff, nor for a deviation from a power 
law at lower energies. The detection of very high energy y- 
rays from Cas A provides evidence of the acceleration of multi- 
TeV particles in S NR shocks and their visibility in gamma-rays 
iDrurvetalJ ( 11994 . 

Significant efforts have been made for the theoretical mod- 
eling of Cas A's multi-frequency emission, including that at the 
highest energies. The effect of an energy-dependent propagation 
of relativistic electrons in a spatially inhomogeneous medium 
has been used in order to interpret the r adio emission from th e 
region and define its electron content (Atova n et al.l (l2000al) ). 
The variations in brightness in the radio band is so complex that 
a multi-zone model was used: distinguishing between compact, 
bright spectrum radio knots and the bright fragmented radio ring 
on one hand, and a diffuse plateau on the other A three-zone 
model with a magnetic field decreasing from its highest value 
in the compact zones putatively related with acceleration sites, 
to a lower value in regions surrounding the shell, to yet a lower 
value in the neighborhood has been found to reproduce the radio 
data, with a magnetic field around and below 1 mG. The fluxes 
at TeV energies, due to Bremsstrahlung and inverse Compton 
radiati on of the sam e relativ istic electrons have also been com- 
puted d Atovan et alj (l2000bl) ) and, albeit the parameters allow a 
large range of possible fluxes, the overall shape of the spectrum 
remains similar, showin g a steep cutoff fo r multi-TeV energies 
(see, e.g.. Figure 7 of Atov an et all (|2000b)). This cutoff is not 
seen in HEGRA and/or MAGIC data, disfavoring a leptonic ori- 
gin of the radiation. Vink and Laming (2003) also studied multi- 
zone models for Cas A, assuming no difference between zones 
other than in their magnetic field. They found that an IC origin 
of high energy fluxes would be possible but only for low val- 
ues (within the range allowed to be consistent with radio and X- 
ray observations, see e.g., Vink and Laming 2003, Hwang et al. 
2004) of the magnetic field and high, far-infrared photon density. 
The generally high values of the magnetic field necessary to ex- 
plain the multi-frequency observations makes it likely that TeV 
emission from Cas A is then dominated by pion decay (Atoyan 
et al. 2000b, Vink a nd Laming 2003). 

iBerezhko et al.l (1200 3) applied a non-linear kinetic model of 
cosmic-ray acceleration to describe Cas A, ignoring the role of 
any small scale inhomogeneities for the production of the very 
high energy particles and considering the whole SNR blast wave 
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Fig. 4: Spectra of Cas A as measured by MAGIC. The shaded area 
around the 0.65 TeV detection shows the la statistical error range un- 
der the assumption of a E"" power law spectrum. The upper limits given 
by Whipple, EGRET and CAT are also indicated, as well as the HEGRA 
detection. The MA GIC and HEGRA sp ectra are shown in the context 
of the model by B erezhko et al.l ( l2003h . Both hadronic (tt" with and 
without an energy cutoff) and leptonic (NB and IC) y-ray emission are 
shown. The normalization of the pion decay spectrum can be taken as a 
free parameter. 



as the main relativistic particle generator. Figure |4] represents 
the expected integral y-ray flux components from non-thermal 
Bremstrahlung NB, IC scattering on the background radiation 
field (cosmic microwave -i- optical/infrared), and hadronic colli- 
sions of CR protons with gas nuclei, respectively, for this model. 
The pion-decay y-ray flux presented in Figure|4]-with and with- 
out an exponential cutoff" at 4 TeV- was calculated with a renor- 
malization factor of 1/6 (i.e. this factor takes into account that 
not all the SNR shock efficiently injects and accelerates cosmic 
rays). This emphasizes that the normalization of nucleonic pre- 
dictions of y-rays is to be considered a free parameter, within 
certain reasonable boundaries. The predicted slope for the dom- 
inating nucleonic-produced y-rays (that dominates, even when 
all possible uncertainties leading to an increase of the leptonic 
emission are included) is hard in the range of interest, as shown 
in Figure |4] perhaps too hard already to provide a good fft to 
the new MAGIC data at low energies. Higher and lower energy 
measurements, and a better signal to noise ratio for the spectrum 
determination of such a weak source, are still needed for a deff- 
nite answer. 
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